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Abstract 

J> ' In this paper we give a central limit theorem for the weighted quadratic vari- 

. ations process of a two-parameter Brownian motion. As an application, we 

Q . show that the discretized quadratic variations X]l=i Sj*^! I'^jj'^l^ of ^ two- 

I parameter diffusion Y = {Y(^s,t)){s,t)^[o,i]^ observed on a regular grid Gn is an 

asymptotically normal estimator of the quadratic variation of y as n goes to 
infinity. 
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1 Introduction 

Many statistical properties of stochastic processes can be deduced from their weighted 
p-power variations. For a one parameter process {Zt)t£io,i] observed at regular times 
{i/n, < i < n}, this quantity is defined as. 
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For example, the study of the power variations process has been used by Barndorff- 
Nielsen and Shephard in [5, 6] to solve some financial econometric problems (applica- 
tion to econometrics are also given in [7]). These theoretical results were also used in 
several fields of application such as the estimation of the integrated volatihty (see for 
example [1] and references therein), testing for jumps of a process observed at discrete 
times like for example in [3]. 

In this paper we give a central limit theorem for the weighted quadratic variations pro- 
cess of a two-parameter Brownian motion. More precisely consider a two-parameter 
Brownian motion W = iW(s,t))(s,t)e[o,i]^ and a deterministic and regular enough func- 
tion / : M ^ M, we show that 

(1.1) 

where S is a two-parameter Brownian motion independent of W and A.ijW denotes 
the increment of the process W on the subset Ajj := ^] x ^] of [0, 1]^ 
defined by 

AijW := + - W(i^ iN - (1.2) 

Vn'ri/ VTi'n/ Vn'ri/ Vn'n/ 

The notation law{S) used above in (1.1) means that the convergence is in the sense 
of stable convergence in law in the two-parameter Skorohod space. Furthermore we 
stress that the limiting process is defined on an extension of the considered probabil- 
ity basis. Note also that usual techniques of proof used in the one-parameter setting 
are no longer suitable to the two-parameter case. For example the Ito formula for 
two-parameter diffusion processes cannot be applied as in the one-parameter setting 
due to the presence of an additional term. Consequently we chose to replace the 
usual stochastic calculus by the Malliavin calculus which is valid in general Gaussian 
context. 



As an application we deduce a central limit theorem for the quadratic variations 
process of a two-parameter diffusion Y = {Y(^s,t)) (s,t) eio,i]'^ observed on a regular grid 
Gn which allows us to construct an asymptotically normal consistent estimator of the 



2 



quadratic variation of Y. More precisely, consider a two-parameter stochastic process 
(^(s,t))(s,t)e[o,i]2 defined by, 

Yis,t) ■■= I (y ^{u,v)) dW(u,v) + / M(„,,) dudv, {s,t) e [0, l]^ (1.3) 

JlO,s]x[0,t] J[0,s]xlO,t] 

where (W^(s,t))(s,t)G[o,i]2 is a two-parameter Brownian motion, a : M — M is a bounded 
sufficiently smooth deterministic function and {M(^s^fj) i^s,t) €[0,1]'^ is a continuous adapted 
process. Assume that {Y(^s,t)){s,t)e[o,i]'^ is observed on the regular grid 

Gn ■= {{i/nj/n), Q < ij < n}, 

let 

[ns] [nt] 

i=i j=i 

and 

Cis,t) := / rfwrf^;, (s,t) e [0, l]l (1.5) 

7[0,s]x[0,t] 

Using (1.1) we show in Lemma 4.4 that. 



["■] [«•] 



Jaui(5) 
n— ►oo 



J2 E i^^'^-^i' - / i^M) ^^^^ (1-6) 

i=l j=i ^[o,-]x[o,.] / 



which is then used to prove that the consistent estimator V" of C is asymptotically 
normal (c/. Proposition 4.2). 

Similar results have been recently established in the one-parameter setting [1, 12, 14, 
15], let us mention some of them. Consider a one-parameter semimartingale {Zt)te[o,i\ 
observed at regular times {i/n, < i < n} with 

Zt^Zo+ [ a {B,) dB,+ [ b {B,) ds, t e [0, 1] (1.7) 
Jo Jo 

where {Bt)t£[o,i] is a standard Brownian motion and cr : M — >■ M and 6 : M — >■ M are 
sufficiently regular deterministic functions. 
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Gradinaru and Nourdin have shown in [12] that, 




(t'^{Bs) ds 



n— >oo 



law 



V2 [a\s,p('^)dpP, 



(1.8) 



where /^^^^ and are two independent Brownian motions. Note that the conver- 
gence obtained in (1.8) hold in the Skorohod space. Property (1.8) has been used by 
Gradinaru and Nourdin in [12] to construct of a goodness-of-fit test for the integrated 
volatihty (their results are even more general since they can be applied to diffusion 
processes of the form (1.7) where the diffusion and the drift terms depend on the 
observed process Z and not only on B). 

In [14, 15], Jacod proved functional limit theorems similar to (1.8) in a larger set- 
ting than in (1.7) since he considered quite general functions of the increments and 
the process {Zt)te[Q,i] was supposed to belong to the class of Ito semimartingales which 
contains some non-continuous processes and Levy processes. We refer to [1] for similar 
results established by Ait Sahalia and Jacod. We also mention that Nourdin in [20] 
and Nourdin, Nualart and Tudor in [21] have studied weighted power variations of 
a one-parameter fractional Brownian motion. Furthermore Nourdin and Peccati in 
[22] have investigated the asymptotic behavior of weighted p-power variations for the 
iterated Brownian motion. 

We proceed as follows. First we recall in Section 2 some elements of stochastic anal- 
ysis of two-parameter processes. Actually we present some definitions concerning 
stochastic calculus of two-parameter processes taken from [13] and the definition of 
the two-parameter Skorohod space initially introduced in [23] and in [31]. Secondly, 
in Section 3 we establish the central limit theorem (Theorem 3.1) for the weighted 
quadratic variations process of the two-paramctcr Brownian motion briefly presented 
in (1.1). As an application we prove in Section 4 that the consistent estimator (1-4) 
of the quadratic variation C (1.5) is asymptotically normal (Proposition 4.2). Finally 
we present in an appendix (Section 5) some background on set-indexed processes, 
extension of probability bases and on the Malliavin calculus for the two-parameter 
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Brownian motion which are used in Sections 3 and 4. 

2 Stochastic analysis of two-parameter processes 

In this section wc recall some definitions of two-parameter stochastic analysis which 
will be used in Sections 3 and 4 and we present the two-parameter Skorohod space 
introduced in [23] and [31]. 

Some elements of two-parameter stochastic calculus 

Let {0,,J^, (•^z)z6[o,i]2,IP) be a filtered probability space. 
We denote the partial order relation :< on [0, 1]^ defined by, 

z' diz^{s' <s and t' < t), z' = {s\t'), z = {s,t). 

We also define the strong past information filtration on {fl, JF, P). 

Definition 2.1. Let z = {s,t) in [0, 1]^. 

s'<s or t'<t 

Until the end of this paper we assume that the following commutation condition hold. 
This property is a conditional independence property (CI in short) and corresponds 
to the condition (F4) of [11]. 

Assumption (CI): 

The filtration (J-'z) z£[o,i]2 is supposed to satisfy the (CI) condition i.e. for all z = (s, t) 
andz' ^ {s',t') in [0,1]^ 

E []E[-|:FJ |JF[o,z]n[o,0']] = E [•|JF(sAs',tAt')] ■ 
Definition 2.2. An {Tz)ze[o,i\'^-O'dapted process {Yz) z£[o,i\'^ is said to be 
i) a martingale if for every z and z' in [0, 1]^ such that z <z' 

W^z'\Tz\ = Yz. 



5 



ii) a strong martingale if for all z and z' in [0, 1]^ such that z < z' 



E[rMi^;] = o, 

where y\z,z'\ denotes the increments ofY on the interval [z,z']. 

As an example, we mention the two-parameter Brownian motion {Wz)zelo,i]'^ is a strong 
martingale with respect to its natural filtration and a centered Gaussian process with 
covariance function, 

nWis,t)W^s',t')] = {sAs'){tAt'), {s,t),{s',t') e [0,l]l 
Skorohod space V{[Q, l]^) 

In the one-parameter setting, Skorohod introduced in [30] four topologies known as 
Ji, J2, Ml and M2. The topology M2 is the weakest of theses topologies in the 
sense that convergence of a sequence (a;„)„ of functions on [0, 1] to x for Ji, J2 or 
Ml consists in the convergence of to x in M2 plus some additional conditions. 

The M2 topology has been extended to the general setting of set-indexed functions by 
Bass and Pykc in [9] whereas the Ji topology has been extended to multiparameter 
functions by Neuhaus and Straf respectively in [23] and [31]. The two-parameter 
Skorohod space (relative to Ji) introduced in [23] and [31] is denoted by ©([0, 1]^) 
and give an equivalent to two-parameter functions of the notion of cadlag functions 
on [0, 1]. The set "^([0, 1]^) can be equipped with a metric d which makes it a Polish 
space and we denote by £2 the Borel u-algebra on {T>{[0, 1]^), d). Note that as in the 
one-parameter setting the Ji topology is stronger than the M2 topology. Furthermore 
compact sets (relative to Ji) on {'D{[0, 1]^), d, £2) can be described thanks to a modulus 
of continuity w which enables us to use techniques described in [10] for one-parameter 
functions. We conclude this section by giving the definition of w. Let / : [0, 1]^ — >■ M 
be an element of I^([0, 1]^) and 5 > we define w{f, S) as, 

w{f,S):^ sup \f(s,t)-f{s',t')\, (2.1) 

\\{s,t)-{s',t')\\<S 

where \\{s,t) - {s' ,t')\\ := max{|s - |t - for {s,t), {s' ,t') G [0,1]^. 
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3 Central limit theorem 



In this section wc state and prove the functional limit theorem (Theorem 3.1) which 
will allow us to show in Section 4 that the consistent estimator V"' (1.4) of the 
quadratic variation C (1.5) is asymptotically normal (Proposition 4.2). 

Let / : M ^ M be a bounded and measurable deterministic function. Let a two- 
parameter Brownian motion W — (VF(s,t))(s,t)e[o,i]2 defined on a probability basis 
B := {n,J^,{T^,^t)){s,t)e[o,i]^,'^)- Let also 

Ci,j:^nf(w^^_^i^yj (^\A,jWf-^y l<i,j<n, n>l. 

The re-normalized weighted quadratic variations process X"- = (-^(^ t))(s,t)e[o,i]2 is 
defined as, 

[ns] [nt] 

^rs,*)^=EE^M' (s,t)e [0,1]'. (3.1) 
i=i j=i 

Stable convergence in law has been introduced by Renyi in [28] and in [29] . It requires 
some particular care, since here the limiting process X is not defined on the probability 
basis B — {Q, T, (J^*)2g[o,i]2, P) on which the X", n>l are defined but an extension 

B:- {Sl,T, (J^.).6[o,i]2,P) of B. 

Theorem 3.1. Assume that the deterministic function / : M ^ M considered above 
is bounded. Then {X"')n>i defined by (3.1) converges !F -stably in law in the Skorohod 
space {T>{[0,l]'),d, £2) to a non-Gaussian continuous process X presented below in 
the proof by (3.2) defined on an extension of the probability basis B. 

Proof. Let us first describe the extension of B on which the limiting process X is 
defined. 

We denote by B' :— {Q! ,T' , (^^)2g[o,i]2, P') the two-parameter Wiener space, that is 
Q! :— C°([0, 1]^) is the space of real-valued continuous functions on [0, 1]' vanishing 
on the set {{s, t) e [0, 1]^, s = or i = 0}. Then P' is the unique measure on {Q', T') 
under which the canonical process (-B0)zg[o,i]2 on defined by, 

B^{J) := J{z), J en', ze [0, 1]', 
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is a standard two-paramctcr Brownian motion. 

Let the extension B := {T*)ze[o,ifj^) defined as, 

(_:F^)2g[o,i]2 := {y\p>zT*p (g) J^^);jg[o,i]2, 
F{dw,dy) := ¥{duj)F'{dy). 

We will denote by E (respectively E) the expectation under P (respectively under P). 
On B we define (-'^z)ze[o,i]2 as, 

X,{u;,u;') := V2 ^ /(W^p(u;)) ciS,^ (u;'), e [0, 1]'- (3.2) 

The process X is a jF-progressive conditional Gaussian martingale with independent 
increments on B, which means that X is an {^z)ze[o,i]'^- adapted process such that for 
P almost cu in fl, X{u!, •) is a Gaussian process on B' with covariance function 

Er [X{si,ti)i^y ■)^{s2,t2)i^^ ■)] 
= 2 /" f {W,) {u) dp, (si, t,), {S2, t2) e [0, 

^[siAS2,SlVS2lx[tiAt2,tlVt2] 

Note that B is clearly a very ^'ood extension of i3 in the sense of Definition 5.4. 

Since (^^([0, 1]^), d, C2) is a Polish space, by [17, Proposition VIII. 5. 33], ^-stable con- 
vergence in law holds if for every random variable Z on (Q, P) the couple {Z, X'^)n 
converges in law. Adapting an argument presented in the proof of [17, Theorem 
VIII. 5. 7 b)], the convergence in law of a such couple (Z, will be obtained as 
follows. First we give a tightness property for the sequence {X'^)n (relative to the Sko- 
rohod space (I^([0, 1]^), c?, £2)) and then we make an "identification of the hmit"wa 
^-stable finite-dimensional convergence in law to X. Recall that the latter prop- 
erty means that for every integer m > 0, for every continuous and bounded function 
-0 : M'""'"^ — >• M and every elements ;2o, . . . , in a dense subset of [0,1]^, 

E [Zi^ (X"(zo), ^"(^1), . . . , X^{zm))] E [Z^jj iXizo),Xizi), X(^^))] . (3.3) 

n— >oo 

The proof is decomposed in two steps. In Step 1) we show that {X"-)^ is tight in 
{'D{[0, 1]^), d, C2) and in Step 2) we prove the ^-stable finite-dimensional convergence 
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in law to X. 



Step 1) 

We show the sequence {X'^)^ is tight in the Skorohod space ('^([0, 1]^), d, £2)- 
A complete description of (I>([0, 1]^), £2) can be found in [23]. In particular it is 
shown in [23] that the set of conditions (3.4) and (3.5) is necessary and sufficient for 
the sequence (X")„ to be tight in (^'([0, 1]^), ci, £2), 



(Xq )n converges in distribution, 
limlimsupP[w(X",5) > e] = 0, £ > 0, 

6 >'0 ri. — vnn 



(3.4) 
(3.5) 



where w is defined in (2.1). Property (3.4) is clear since for every n>l X^ — Xq — 
0, P-a.s.. We will show (3.5) using a method from [10, p. 89]. 

Let e > 0, 5 > and n > 1. Let m := [|] and v := [^] . We consider on [0, 1]^ the 
rectangles R,, := f ] x ^] , E {!,■■■ where rm := im, 1 < 

i < V and = n. With this notation the length of the shortest side of the rectangles 
Rij is greater than 6 and v < 2/6. We can adapt the proof of [10, Theorem 7.4] to 
our case and we have. 



P[w(X",5) > 3£] < ^ 

i=l j=l 



sup 



X': - xr 



n ' n , 



> e 



(3.6) 



Let us give some notations. For (A;, j) e {1, . . . , n}^ let 



k I 



1=1 j=i 



that is XVj^i^ = nSk^i- For z in Ri j we write S^^i Sk,i — S mi_i mj_i . Using these 
notations we can write (3.6) as. 



P[w(X",(5) > 3£] < EE^ 

i=i j=i 



sup 

mi-i<k<mi, mj-i<l<mj 



> 



n 



We will now use [10, Section 10] which provides maximal inequalities for partial sums 
of non-independent and non-stationary random variables. For i,j fixed as above, we 



9 



re-index the random variables appearing in S^^-'^ to obtain, 

r]{ij,k,l) 

p=i 

with Tp equal to some divided by n and r){i,j, k, I) is an integer. 

Let two integers a < (5. Since / is supposed to be bounded by a non-random function, 

let it! := sup^jgg non-random. Let K, K denote non-random constants. 



P 



p=a+l 



> A 





p=a+l 



h E n\rp\ 



< 



A^ 

p=a+l 



(3.7) 



(3.8) 



Using [10, Theorem 10.2] and (3.7) we obtain 



P 



max Sj^] > X 

i<k<m,j<l<m ' 



< 



KBSr&_ 



(3.9) 



Now injecting inequality (3.9) in (3.6) we have. 



P[w(X",5) > 3£] < 



< —. — —, since v < 2/0, 



< 



-n 



4(p-i)^4p-2^ Since m = [n5]. 



which leads to (3.5). 



Step 2) 

Here we choose to consider processes X" and X as set-indexed processes and we use 
all the notations and definitions of Subsection 5.1. Consequently the jF-stable finite- 
dimensional convergence in law property (3.3) can be rewritten as follows: for every 
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continuous and bounded function ip, for every elements Cq, . . . , in a dense subset 
of A (see Subsection 5.1 for definitions and notations) and for every random variable 
Z on {fl,J^,F), 

E [Z^lJ {X-{Co),X^{C,), . . . , X-{CJ)] E [Z^|J {X{Co), X{Ci), . . . , X{C^))] . 

n— *oo 

(3.10) 

To obtain (3.10) we adapt [13, Proposition 7.3.7] which allows us to replace ^-stable 
finite-dimensional convergence in law with T-stable semi-functional convergence in 
law that is, for every simple flow (p (see Definition 5.1) the sequence of one-parameter 
processes {X"' o (p)^ converges ^-stably in law to the one-parameter process X o (p. 
Let us make precise this argument. 

Assume that stable semi- functional convergence in law holds. We aim at showing 
(3.10). As in [13, Proposition 7.3.7] since for every n > 1, is an additive process 
(see Definition 5.2) it is enough to prove (3.10) for elements Cq, ■ ■ ■ , Cm such that there 
exists a simple flow ip such that for every i e {1, • • ■ , m}, Ci = ip{i/m) — ip{{i — l)/m). 
Since the sequence of one-parameter cadlag processes (X"o(^)„ is supposed to converge 
^-stably in law to X o(p, and since one can choose a continuous version of X o(p, then 
the projection 7r(o,i/m,...,i) '■ ^([0, 1]) W"'^^ is continuous and by mapping Theorem, 

E [Zi; ((X" o (^)(0), (X" o ip){l/m), ■■■,{X-o (^)(1))] 

Z(u;)'ijj ((X o ip)(u;, x)(0), ■ ■ ■ , (X o ip)(u;, x)(l)) r(dx)¥{du;) 

= E [Z^ {{X o (^)(0), {X o ^)(l/m), • • • , (X o ^)(1))], (3.11) 

where Z and are like in (3.10). 
Consequently relation (3.10) holds since 

X'^id) = (X" o ip){i/m) - (a:" o ip){{i - l)/m), 1 < i < m. 

Using the argument presented above we will now prove ^-stable semi-functional con- 
vergence in law to establish J^-stable finite-dimensional convergence in law. 

Let </? be a simple fiow (we write (p as (p = {pi, P2))- We have to show that the sequence 
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of one-parameter cadlag processes (X" o converges jF-stably in law to the one- 
parameter process X o(p. We give some precisions about the extension of probabihty 
basis we use. We set := {J^^^t))te[o,i],^) and B'^ := {J^^^^^)te[o,i]X)- 

Prom Bip and B'^ we define the probabihty basis B^p :— (.F^(t))tg[o,i], P), with, 

(■^¥'(t))te[o,i] — {^s>t^ip{s) ® •^<^(5))te[o,i]- 

Let n > 1, by Lemma 5.3 and Lemma 3.2 the one-parameter processes ((X"o(p)j)jg[o,i] 
are martingales on the probability basis B^ := {fl,J-', {J-'f''^)te[o,i],^) where, 

•^'"^ ^([„<^i(t)]n-i,[n<^2W]n-i), t G [0, 1]. 

We define also, 

/ •^(Wl(t)]«-^[«'P2(t)]n-l)' ^ ^ [0' 1]' 

which lead to the following probability basis, 

r B'^ := in',r,irrUo,i],r), 

\ Bl-.^ (Q,J^,(J^'%e[o,i],P). 

We will now apply [17, Theorem IX.7.3] (or its triangular array formulation [17, 
Theorem IX. 7. 28]). This result gives conditions insuring ^-stable finite-dimensional 
convergence in law for a sequence of one-parameters martingales to a continuous 
conditional martingale with independent increments by identifying the caracteristics 
of these martingales plus some additional conditions. This identification is realized in 
Lemma 3.3 in which convergences (3.12) and (3.13) can be thought as identification 
of the characteristics whereas properties (3.14) and (3.15) ensure the ^-stable feature 
of the convergence. Consequently from [17, Theorem IX.7.3] and Lemma 3.3, the 
sequence (X" o y?)^ of one-parameter martingales on B^j converges jT-stably in law to 
X o if on the extension B^ of B^ which ends the proof. Note that B^ is a very good 
extension of B^ since i3 is a very good extension of B. □ 

Before turning to estimation results in Section 4 we state and prove Lemma 3.2 and 
Lemma 3.3 which were used in the proof of Theorem 3.1. 
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Lemma 3.2. We use notations of Theorem 3.1 and of its proof. X"- is a strong 
martingale. 

Proof. Let y and x in [0, 1]^ such that y = (1/1,1/2) ^ x = {xi,X2). Let also, 



-pn,* 



(n i[ni/i],n ^[nj/2])" 



We have 



E[x"([y,x])|^;'*] 

[nzi] [na::2] p . 

j=[nj/i] j=[nj/2] 

[nxi] [na;2] 



pn,* 



0. 



□ 



Lemma 3.3. We use notations of Theorem 3.1 and of its proof and in particular 
we denote a flow (p as (p — {(fii,(fi2). For every n > 1, X"' o (f is a one-parameter 
martingale with modified second caracteristics (0, Cx"o<p, ^x^oip) on such that 

i^xno^([0, t] X {|x| > e}) -^0, e [0, 1], £ > 0, (3.12) 

n— >oo 

Cx^o^it) ^ 2 / f {Wp) dp, yt e [0, 1]. (3.13) 

7[0,v.i(t)]x[0,^2{t)] 

Furthermore for every bounded martingale N orthogonal toWoip, 

(X" o 7V)t 0, Vte[0,l]. (3.14) 

n— >oo 

We also have, 

(X'^ o o (^)i 0, Vte[0,l]. (3.15) 

n— ►00 

Proof. Let n > 1. 
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Proof of (3.12): 

In the following /ix^oip denote the jump measure of X" o ip. For lu & fl, 



fJ'X^o<f{uJ, dt, dx) - ^ [^{<^-i(ifcn)}n(i^2 i({/n, l<i<n}))=,i?i.'=."(a;) i^^-: '^^) 

k=l 

+ ^{ip-\kn)}n{v-^{{ln, l<l<n}))'',H2,k,n(^^){dt, dx)] 

n iL 

+ X] X] ^>p^''{k/n),v:^\i/n),H^^'',i,u^^^{dt, dt, dx), 



k=l 1=1 



We denote by i^x"o(^ the compensator of the measure nx^oif- Let ^4 a Borel set in R, 
we have 



[nipi{t)] 

E [l/fi,fc,n6A|jF(j^_i)/„,([„i^2(t)]-l)/n] 

fc=l 

[n(p2(t)] 

+ E [l^f2.i,ng^|J^(([„^i(t)]_l)/„,(i_l)/„)] 

[n<^i(t)] [n<j52(t)] 

+ XI X] ^ [li?3,fc,i,ngA|-^([n<^i(t)]-l)/n,([n<^2(t)]-l)/n]- 
fe=l 1=1 

Let £ > and k,l,n,t as above. Since / is assumed to be bounded by a non-random 
constant let R :— sup^jgj^ 1/(^)1- Denote by C and C some constants. 

4 



1 



[nip2{t)] 



k-l [ri' 



[nip2(t)]-l ^ 
n ' n J 



(*) 1 



n ' n 



g 4 [n>f2it)] 



< 



< 



R^C 



J2 E \{\A,,W\'-l/n' 
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Before giving details about the inequality (*) we deduce of the preceding inequalities 
that, 



P 

< P 



I ^3,k,l,n 
|^l,fe,n| 



> e 

e 

>3 



/ [nyi(f)]-l [ny2(t)]-l \ 
Y n ' n y 



/ [ny>i («)]-! [ny2W]-l N 
Y n ' n y 



+ P 



2J,rt| 



>3 



Y n ' n y 



+ P 



Which leads to (3.12). Now we give some details about inequality (*). 
Let 1 < h, h, h < h < [n(p2{t)] — 1- We have 



IE 
E 



[n<p2(t)]-l \ 
\ n ' 71 J 

■ ■ -^^4/ E (|Afc,,,iy|^ - 1/n^) 



\^ n ' n y 



= nE 
= 0. 



■J' ( fc-i 



[wy2 (*)]-! 
n ' n 



Proo/ 0/ (3.13): 

Cxno<p{t) — {X'^ o X" o is the compensator of [X" o X" o with respect to 
{J^''^)te[o,i] (see for example [17, Proof of Proposition II. 2. 17 b)]) and we have (c/. 
[17, (1.4.53)]), 

[n<pi{t)] [nip2{t)] 

[X"o(^,X"ovp]t= J](X"o(^)(s)-(X"o<^)(s_)= J] 5^4, te[0,l]. 

0<s<t i=l j=l 



Consequently, 

Cx"o:f{t) 



[n<pi{t)] 

E [(i7^''''")^|j^(fc-l)/n,([nvP2(t)]-l)/n] 



fc=l 
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In<fi2it)] 

+ ^ E [{H^'^'''y\J^^(^lri^^(t)]-l)/n,(l-l)/n)] 
1=1 

ln<fii(t)] [nip2{t)] 

+ XI ^ [(^^''''''")^l-^(["¥'i(t)]-l)/n,([nV2(t)]-l)/n]. 

k=l 1=1 



We can show this sum is equal to 

[nipi{t)] [nip2{t)] 



(3.16) 



i=l 3=1 



since terms of the form E 



the form E 



kd 



y n ' n y 



vanish for j < I < [n(f2{t)] using 



the same type an argument described in the proof of (3.12). Furthermore terms of 



Y n ' n y 

We deduce (3.13) from (3.16). 



are given by, 



E 



E 



\Al^W\'-l/n' 



Proof of (3.14): 

Let iV be a martingale orthogonal to W oip. Without loss of generality we can assume 
there exists a strong martingale on i3 orthogonal to W such that NW is a strong 
martingale and such that N = N o ip. Let n > 1 and t in [0, 1]. We have, 



i=i j=i 

[nipx(t)\ [nip2{t)] 



= E E/'"W. 



i=l j=l 



V n ' n / 



E 



A,,,-iy|' - l/n^) A,,,iV|^.^^ 



Let 1 < 1,3 < n. We use a technique presented in [15, Lemma 6.8]. For z in [0, 1]^ 
we define Uz := E [|Ajjiy|^ — l/n^lJF^]. {Uz)zy{{i-i)/nXj-i)/n) is a martingale for the 
filtration generated by — U/i-i Using the representation |Ajjl4^|^ — ^ = 

I2 ^l^i^) ^ ^ multiple stochastic integral (see for example [26, Section 1.1.2] or [19]) 



16 



we have by [26, Lemma 1.2.5] that for zy{{t- l)/n, (j - l)/n), = h .l[o,z]) ■ 
Prom [18] there exists a adapted process ($(s,t))(s,t)e[o,i]2 such that 



U(s,t) = U^i^J^^ + 



'[{i-l)/n,s]x[{j-l)/n,t] 

Define N'^ = — Nti-i z ^ {{i — l)/n, (j — l)/n). This process is orthogonal 
to {Uz)z>{{:i-i)/n.{j-i)/n)- Consequently using a characterization of orthogonal two- 
parameter martingales given in [11, Proposition 1.6] we have that 



E 



3z 

n ' n 



0. 



A straightforward computation gives that. 



E 



\ n ' n / , 



E 



\ n ^ n J . 



Proof of (3.15): 

Let n > 1 and t e [0, 1]. We have, 

[rvpx(t)\ [nip2(t)\ 

i=i i=i 

[riipi(t)] [riip2(t)] 

= 0. 



71 ' n 



□ 



4 Estimation of the quadratic variation and asymp- 
totic normality of the estimator 

In this section we prove an asymptotic normality property (Proposition 4.2) for the 
consistent estimator V" (see (1.4)) of the quadratic variation C (defined in (1.5)). 



Consider the following two-parameter stochastic process. 



Yis,t) ■■= I (T{Wp) dW, + f Mp dp, {s, t) e [0, 1]^ 

J[0,s]xlO,t] ^[0,s]x[0,t] 



(4.1) 



17 



defined on a probability basis (■^(s,t))(s,t)e[o,i]2) P) • Until the end of this paper 

we assume that ()g[o^i]2 is a continuous and (^(s,t))(s,t)e[o,i]2-adapted process. 

Let the following assumptions which will used in the results presented below. 



Assumption (Rl): 

The function a{-) and its derivatives are assumed to be bounded by non-random con- 
stants and cr(-) is assumed to be at least in C^. 

Assumption (R2): 

The function a(-) and its derivatives are assumed to be bounded by non-random con- 
stants and a(-) is assumed to be at least in C^. 



Let us define the quantities we will study. 
Let for n > 1, 



[ns] [nt] 



— n 



J2 \\,jyf -I {Wiu,v)) dudv . (4.2) 

i=l j=l J[0,s]x[0,t] I 



Proposition 4.1. The estimator V"' defined in (1-4) of the quadratic variation C 
(1.5) is consistent that is for every {s,t) in [0, 1]^, 



^ ' ' n—*oo 



Proof. Fix (s,t) in [0,1]^. We show (V^"j))n converges to C(^s,t) in LF'iyi^J^,. 
enough to show that 

[ns] [nt] . . 

tends to zero in L^(0, J^, P) as n goes to infinity. We have, 



It is 



E 



[ns] [nt] 

i=\ j=i 
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[ns] [nt] 



i=i j=i 



n n , 



E 



n 



- 1 



where R — supj.g]g 

We state and prove that the estimator of C is asymptotically normal. 



□ 



Proposition 4.2 (Asymptotic normality). Let for {s,t) in [0,1] 5'"^^^ be 

[ns] [nt] 



[ns] [nt] 

i=i j=i 



n^QO y o 



Under assumption (R2) for {s,t) fixed in (0, 1]^ we have, 

Proof. Using a localization argument, only the stochastic integral part of Y gives a 
contribution to the limit so we assume M = in (4.1). The main argument of the 
proof is the convergence in law of (5*", y")„ to (5", X) where 5* is defined by, 

Sis,t) := 3 / {W,) dp. 

JlO,s]x[0,t] 

Actually assume this convergence hold. Since (x, y) i— > x'^y is continuous on x M 
we have for every (s, t) in (0, 1]^, 



law^ 
n— >oo V 3 



2 I[o,s]x[o,t]^'WdBp 



(/mxM^MW^p) dp 



1 • 

2 



Computing the characteristic function with respect to the probability measure P we 
can show that 

2 V 3 



(/mxIM^'TO dp) 
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We have now to show that {S"',Y"-)n converges in law to {S,X). The key point is 
the JF-stable convergence in law of {Y"-)n to X obtained in Lemma 4.4 stated and 
proved at the end of this section. Using a result of Aldous and Eagleson (presented 
in [4]) concerning stable convergence in law if (5'")„ converges in P-probability to 
S then (S"", Y'")n converges in law to {S, X) (and the convergence is even ^-stable 
convergence in law). 

Let us finally show that (5'")n converges in L^(n,.F, P) to S. 
First we show that 

2n 



E 



n 



i=i j=i 

Actually for every {s,t) in [0, 1]^, 

[ns] [nt] 



[ns] [nt] . 



0. 



E 



i=\ i=l 

[lit] 



i-1 



"'EEe 

i=\ j=i 



i-l j 

n ' n / / 



c 



where C is a constant. 

Using a Riemann approximation for integrals we have that 



EE ^y~^^ 



=1 i=i 



3 / cj^{W,)dp, V(s,t)e[0,lf. 

'[0,s]x[0,t] 



The proof is finished if we can show the estimate (4.3) 



E 



n 



1 j=i 



0. 



(4.3) 



Using the same techniques (successive Malliavin integrations by parts, estimates of 
the form (4.9)) and (4.10) we obtain (4.3). 



□ 
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Using classical techniques the asymptotic normality property of V" enable construc- 
tion of confidence interval for C(s,t) for every {s, t) in [0, 1]^. 

Remark 4.3. It seems that conditions on the regularity on cr(-) imposed in Theorem 
3.1 and in Proposition 4-2 are too strong. These type of regularities are imposed by 
the Malliavin calculus techniques we use. However different techniques used in the 
one-parameter case like for example the ltd formula are much more difficult to use 
in the two-parameter setting. In contradistinction, the Malliavin calculus presents a 
large range of properties valid in general Gaussian context. 

In order to prove the asymptotic normality of V"' in Proposition 4.2 we have used the 
following lemma. 

Lemma 4.4. Let X be the process defined in (3.2) where f in (3.2) is replaced by a^. 
Under Assumption (Rl), (V'")n>i defined by (4-2) converges -stably in law in the 

Skorohod space (V{[0,l]'^),d, C2) to the non-Gaussian continuous process X defined 
on the extension B described in the proof of Theorem 3.1. 

Proof. Using a localization argument the finite variation part of Y has no contribution 
in the limit. So we can assume that M = 0. 

From Theorem 3.1 with f = a^, the process converges JF-stably in law to X 

with, 



To conclude the proof we show that is equal to X" plus a term r„ which become 
negligible when n to infinity. More precisely using the notations 





r„ can be decomposed as, 



(4.4) 
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where 



ri? (s, t):=- jl^y^ ^2 iWp) dp. 

Using a standard argument of the form [10, Theorem 3.1] as converges JF- 

stably in law to X it is enough to prove that wsup(5_^)g[o,i]2 kn(s,t)| converges in 
probability to zero to obtain the J^-stable convergence in law of {Y"')n to X. We use 
the decomposition (4.4) and we show 

n sup |rn(s,t)^^^| ^ 0, (4.5) 
(s,t)e[o,i]2 

n sup i)(^)|^0. (4.6) 

(s,t)e[o,i]2 

Proof of (4.5): 

Using Burkholder's inequality for two-parameter martingales (see Remark 2 of [24]) 
it is enough to show that 

E 



Vi,j 



•3 

where C is a constant. 

The tool used here is the MaUiavin calculus (see Appendix 5.3) and especially the 
Malliavin integration by parts formula (5.3). The main problem comes from the 
computation of E [|?7ij|^] . First we express r/j j as, 

j a {Wp) - a (W^^,^^^ dW^ (^j a{Wp) + a (w^i^i^-^^ dWp 
= 5{u)5{v), 
where 

S{u) denotes the Skorohod integral of the process {u(^s,t)){s,t)elo,i]'^ which coincides since 
u is adapted with the Ito stochastic integral. Consequently, 

= mu) {6{u) 6{vr)] 

= E[{u,D{S{u)5{vf))L^^lo,ir)], by (5.3) 



{s,t) e [0,l]l 
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[ lE[u^s,t)D^s,t){^{u)S{vf)]dsdt 

JAij 

]E[u^,,t)5{vfD^sd^{u))] dsdt 



(4.8) 



+ 2 E S{u) 6{v) D^.^Siv))] dsdt, 



where the last equahty is deduced from (5.4). We can then apply successive Malliavin 
integration by parts to each term of the right hand of (4.8). Using the "Heisenberg 
commutativity relationship" (5.5) and estimates (4.9) and (4.10), we obtain 

C 



rr 



where C denotes a constant different from the one presented in (4.7). 



E 



\ n ^ n / 



pe\j, A; = 2,4,6,i 



(4.9) 



where the denotes some constants. 

Let (/? : R — > R a deterministic function in C^, 



= if' (W(,-l 3-1 \ 
\ \ n ' n J 

+ «'"('*'(-.^))('^'-'*'(-.^)) 



Wp - Wf^^ 

\ n '' n / 



(4.10) 



where f" is a negligible term. 

Finally note that the constant C in (4.5) depends on the following deterministic 
bounds 

suPpe[o,i]2 sup(,^i)g[o,i]2 \Dpa {W^s,t)) \ = sup^-eR 
SUP(p,a)e([0,l]2)2 SUP(3_i)g[o,l]2 \De,Dpa {W^s,t)) I = sup^^M 

where D denote the Malliavin derivative (see (5.2)). 
Proof of (4.6): 

This convergence result is obtained from Markov inequality and by standard argu- 
ments in numerical analysis. □ 
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5 Appendix 



In this section wc present some definitions and results used in Sections 3 and 4. First 
we provide some background on set-indexed processes and on extensions of probability 
bases. Finally we briefiy present the Malliavin calculus for two-parameter Brownian 
motion, including the Malliavin integration by parts formula which has used to obtain 
the estimates of the previous sections. 

5.1 Set-indexed processes 

In the following definition it will be convenient to think of two-parameter processes 
on [0, 1]^ as set-indexed processes on A := {[0,2;], z G [0, 1]^} where Y[o,^] := and 
^[o,z] '■= ^z, z in [0, 1]^- We will use indifferently one of these two points of view. 

Definition 5.1. ([13, Definition 7.3.1]) Let A{u) he the set of all finite unions of 
sets from A. A simple flow </? is an application </? : [0, 1] A{u) which satisfies the 
following properties, 

i) (p is increasing, 

a) ip is continuous, 

III) <^(0) = 0, 

iv) for some k in N there exists some increasing functions (pi : [0,1] ^ A , i — 
1, . . . ,k such that for ^ <s<|; 1 < i < k 



Definition 5.2. ([13, Definition 1.4-5]) Let C denote the set of elements of the form 
A\B with A in A and B in A{u) (the set of all finite unions of elements of A). 
A process {Xz)zeio,i]2 identified with (Xa)ag^ is called additive if for every C,Ci,C2 
in C with C = C^U C2 and Ci H C2 = we have. 




Xc — + Xc2 ■ 
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We recall a particular case of [13, Lemma 5.1.2]. 

Lemma 5.3. Let X be a strong martingale and let ip he a A{u)-valued simple flow 
defined on [0,1]. Then the one-parameter process {X^p|^s)) s&[o,i] is a martingale with 
respect to the filtration {J-'^(^s))sgIo,i]- 

5.2 Extension of a probability basis 

We give the definition of a very good extension of a probability basis which has been 
introduced in [16]. The following definition is an adaptation of [17, Definition II. 7.1]. 

Definition 5.4. A prohaMlity basis B = ^, {^z)zeio,i]^T^) is an extension of the 
probability basis B — {fl,J-, {J-z)ze[o,i]'^T^) if there exists an auxiliary probability basis 
B' = {Vt\P, (.F^);,e[o,i]2,r) such that 

jr = jr jT', 

= r^p^z^Fp ®r„ ze [0, i]2, 

P(du;, du;') = ¥{du;)Q^{du;') , 

where Q^{duj') is a transition probability from into {Q',J^'). 

This extension is called very good is for every z in [0, 1]^ and for all element A' in T'^ 

Lo 1-^ Qa;(^') is cqual F-a.s. to an Tz-measurable random variable. 

Notation: 

// {Xz)zelo,i]'^ is a stochastic process on B we will denote by {Xz)z£io,i]'2 again the 
stochastic process defined on B by 

Xz{cu,u') ■.= Xz{uj), z E [0,lf, {to, to') e n. 

Lemma 5.5. As in [17, Lemma n.7.3], under Assumption (CI) an extension B 
is very good if and only if every martingale on B is a martingale on B. 

The proof of this Lemma is similar to its one-parameter counterpart [17, Lemma 
II.7.3]. 
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5.3 Malliavin calculus for two-parameter Brownian motion 

Let {W(^s,t)){s,t)€[o,ip be a two-parameter Brownian motion defined on a probability 
basis (fl,J-', {^z)z(z[oi]2 ,Pj. This process is a centered Gaussian process whose co- 
variance given by, 

E [W^s,t)W^s',t')] = A s') {t A f), {s, t) e [0, l]^ is', t') e [0, 1]^. 

The MaUiavin calculus for general Gaussian processes has been described in [26] and 
the reader can refer to it for a complete explanation about this topic. Here we give the 
definition of the Malliavin derivative and we present the integration by parts formula 
which is hardly used in Section 3. 

Definition 5.6. Let S be the space of random variable F of the form 

F = fiWih,),...,W{hr,)), (5.1) 

where hi is an element of {[0,1]^, dz) and W{hi) denotes the stochastic integral 
W{hi) := /[Q hi{z) dz for i = 1, ■ ■ ■ ,n and / : M'* — > R is infinitely continuously 

differentiable. 

For F of the form (5.1) we define the Malliavin derivative DF of F as the following 
L^([0, 1]^, dz) -valued random variable, 

n 

DF:=J2dif{W{hi))hi. (5.2) 
1=1 

Here we give the Malliavin integration by parts formula ([26, Lemma 1.2.1]). 

Lemma 5.7. Let F in S and {u(^s,t)){s,t)e[o,iP domain of the divergence operator 

S. We have, 

E [FS{u)] = E [{DF, 5{u))L2iio,i]2,dz)] ■ (5.3) 

S is called the divergence operator ( or the Skorohod integral of u) and it extends the 
ltd integral since when u is adapted to the filtration generated by the two-parameter 
Brownian motion W , 
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Note also that the Malhavin derivative D is a closable operator (see [26, Proposition 
1.2.1]) and we denote by Dom{D) its domain. Furthermore D satisfy a chain rule 
property that is, for every F and G elements of Dom{D) such that F G belongs to 
Dom{D) we have 

D{FG) ^ F DG + DFG. (5.4) 

We end this section with the "Heisenberg commutativity relationship" which enables 
the computation of the gradient of a Ito stochastic integral, more precisely we have 
for a process u such that the right hand side of (5.5) is well-defined (more details 
about assumption on u can be found in [26, (1.46)]), 

D^s,t)S{u) = u^s,t) + S{D^s,t)u), (s, t) e [0, 1]2. (5.5) 
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